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Abstract 

The chiral phase transition is studied in an extended Nambu-Jona-Lasinio model with eight- 
quark interactions. Equations for scalar and vector quark densities, derived in the mean field 
approximation, are nonlinear and mutually coupled. The scalar-type nonlinear term hastens the 
restoration of chiral symmetry, while the scalar- vector mixing term makes the transition sharper. 
The scalar-type nonlinear term shifts the critical endpoint toward the values predicted by lattice 
QCD simulations and the QCD-like theory. 
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Qualitative properties in quantum chromodynamics (QCD) at high temperatures and 
densities attract much attention. One of the most important recent findings is strong cor- 
relations in the quark gluon plasma just above the critical temperature; it is realized as the 
near perfect fluidity [l|, and the mesonic correlations 14, 51. 

With the aid of the progress in computer power, lattice QCD simulations have become 
feasible for thermal systems at zero or s m all d e mity S0.A thighdemity ,Wve r , lattice 
QCD is still not feasible due to the sign problem. As an approach complementary to the 
first-principle lattice simulations, we can consider several effective models. One of them is 
the Nambu-Jona-Lasinio (NJL) model. Since it was proposed although this is a model 
of chiral symm etry that does not possess a confinement mechanism, this model has been 

A), for example, for analyses of the 

id. 



widely used |9j, |lCj in the mean field approximation Ch '. 



14. 



15, 



critical endpoint of chiral transition ^ , 

Although the NJL model is recognized as an useful method for understanding the chiral 
symmetry breakin g, o nly a few studies were done so far on roles of higher-order multi- 
q uar k interactions Q, except for the case of the si*- q uar k interaction coming from the 
t Hooft determinant interaction [19]. The NJL model is an effective theory of QCD, so 
there is no reason, in principle, why higher-order multi-quark interactions are excluded. In 
the nonperturbative renormalization group method j^j], such higher-order interactions are 
produced as a result of quantum effects in the high momentum region. Such effects should be 
included in the low-energy effective action from the beginning. Thus, it is quite meaningful 
to study effects of higher-order interactions on the chiral phase transition. 

In this paper, we consider an extended NJL model that newly includes eight-quark in- 
teractions and analyze roles of such higher-order interactions on the chiral phase transition. 
It is well known that the original NJL model predicts a critical endpoint to appear at a 
lower temperature (T) and a higher chemical potential (//) than lattice QCD pj and the 



QCD-like theory \2l\ . |22j do. We will show in this paper that a scalar-type eight-quark in- 



teraction shifts the critical endpoint toward values predicted by the lattice simulations and 
the QCD-like theory. 

As for the repulsive vector- type four-quark interaction {qj^q) 2 , it is well-known that it 
makes the chiral phase transition weaker in the low T and high /i region, so there is a 
possibility that the transition becomes a crossover in the region when the interaction is 
strong enough Q|. In this point of view, an absence of the vector- type four-quark 
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interaction may be preferable in the high density region; the absence is also supported by 



works of Refs. 



and 



24 1 . 



On the contrary, in the relativistic meson- nucleon theory |25( , the repulsive force mediated 
by vector mesons is necessary to account for the saturation property of nuclear matter. 
Using the auxiliary field method, we can convert quark-quark interactions to meson-quark 



interactions; for example, see Refs 







27,l2a| and references therein. It is then natural to 



think that there exists a relation between the meson-nucleon interaction and the quark-quark 
interaction in the NJL model, since a nucleon is composed of three constituent quarks each 
of which has a large effective mass as a result of the spontaneous symmetry breaking (SSB) 
of chiral symmetry. In this point of view, the vector-type four-quark interaction (q^^q) 2 is 
indispensable around fi ~ 308 MeV corresponding to the normal nuclear density region. 
Thus, it is expected that the vector-type interaction is sizable in the normal density region 
but suppressed in the higher density region. 

In the relativistic meson-nucleon theory, it is known that nonlinear meson-nucleon inter- 
actions suppress the effective coupling between mesons and nucleons in the higher density 
region [2ia, |2iJ |3iJ, |3J| . It is then strongly expected that a similar situation takes place in the 
NJL model as soon as higher-order multi-quark interactions are introduced. This is shown 
in this paper. 

This is the first work to study roles of eight-quark interactions on phase transitions, so 
we concentrate our analysis on the chiral phase transition and use a simple model with two 
flavor quarks. Effects of the higher-order interactions on the color superconductivity will be 
discussed in a forthcoming paper. 

We start with the following chiral-invariant Lagrangian density with two flavor quarks 

£ = q(ifi- m )q+ g 2 ,o({qq) 2 + (qi^rq) 2 ^ +# 4 ,o {[qqf + (nrq) 2 ^ 

-goAvYq) 2 - 92,2 [{qqf + (^7s^) 2 ) {q^qf + •••], (1) 

where q, mo and g^j stand for a quark field, a bare quark mass and coupling constants. Here, 
we consider only four- and eight-quark interactions by ignoring the higher-order interactions 
denoted by dots in Eq. (JIJ. Furthermore, we disregard interactions including isovector- 
vector current not important in symmetric quark matter, and also does the vector-type 
eight-quark interaction (qj^q) 4 , because the expectation value of the vector current qj q is 
smaller than that of the scalar one qq, unless the chemical potential is quite large. The mean 



field approximation reduces the Lagrangian to 



-MFA 



q{i 
+ 



m )q 



2# 2 ,o + 4^4,o ((qq) 2 + (qh^rq) 2 ^ - 2g 2)2 (q^q) 2 ((qq)qq + (qilbTq)qi^ b Tq^ 
2#o,2 + 2^,2 ((qq) 2 + (qhnrq) 2 ^ (qi^q)qi^q 
-92,o {(qq) 2 + (qil$rq) 2 ) - 3# 4 ,o {(qq) 2 + (qilbTq) 2 ) 2 
+3(7 2 , 2 {(qq) 2 + (qi^qf) {q-fqf + Q^q^q) 2 - (2) 

Below, just for simplicity of the notation, we will omit terms including the pseudo-scalar 
current qij^fq and the spatial components q^iq (i = 1,2,3) of the vector current, since 
their expectation values vanish. It is convenient to introduce two auxiliary mean fields as 



a=(qq), u = (q^ q) . 

Using these auxiliary fields, one can rewrite £mfa as 

£mfa = q[i&- (mo + S B ) + Klo]q - U , 



(3) 



(4) 



where 

E s = - (2g 2fi v + 45f 4i0 cT 3 - 2g 2j2 au 2 ^j , E Y = - (2g 0}2 fv + 2g 2t2 a 2 u^ , 
U = g 2 , cr 2 + 3g 4 ,ocr 4 - 3g 2t2 <7 2 u 2 - g Q)2 uj 2 . (5) 

The thermodynamical potential Q of the system with finite T and /i is then obtained by 



Q(T,ti) = -2N f N c V J 
+VU, 



d 3 p 

(2tt) 3 



E p + - hog (1 + e -0(*W)) + log (1 + e ~^ + ^ 

fj 



(6) 



where (3 = 1/T, /i* = /x + E v , E p = a/p 2 + M 2 and M stands for the effective quark mass 
defined as M = m Q + E s . The corresponding scalar and vector quark densities, p s and p v , 
are also given by 

d 3 p M ( \ oAr Ar f d 3 p 



Ps = 2N f N c j jL^M-fa + n^-i), Pv = 2iV / iV c J -^(n q -n q ) , 



(7) 



where n q = 1/[1 + exp {(3(E P — //*)}] and 



Tin 



[1 + exp {p(E p + ji*)}]. The physical solu- 



tions of a and u satisfy the stationary condition 



Here we have defined four effective couplings, G* a , G* sw and G*^, as |28l 123. 130L 131 ] 



G*„ = = 2g 2fi + 12g 4fi a 2 - 2g 2 ^\ G* va = = Ag 2 , 2 au. (9) 



BE BE 
GL = —q£ = -*92#™ = Gl w = = 2g 0>2 + 2g 2 . 2 o\ (10) 

When det((?*) = G% a G% u + G* a 2 ^ 0, which is satisfied in our analyses below, the matrix Q* 
has its inverse, and then the stationary condition (jBJ) leads to a = p s and u = p v , showing 
the consistency with Eq. (jSJ) [30]. However, the solutions a and uj to the equations do 
not necessarily yield a minimum of Q. The solution a is then determined by searching for 
minima of f2(a,u(a)) in which uj is eliminated with u = p v . 

Identifying qijsrq with the pion field, we can define the effective coupling G* v between 
pion and quark as 



5 3 

lGl ^ = 5q6(^ 5 T q )6q CMFA 



= i (2#2,o + 4#4,oO- 2 - 2fi-2,2^ 2 ) 75-r- (H) 

(qij 5 rq)=0 



In the random phase approximation (RPA), the pion mass M n at T = p = is determined 
with this effective coupling as 



G^itf/filf.M,)' v V (2ir) 3 ^ p 2 + I 2( 4(p 2 + _ y 2 

Similarly, the a-meson mass M a at T = p = is determined as 

. r9 4m + 32o 4 (T 3 , , , 9 , , 

M = — — h 4M 2 . (13) 

This equation indicates that M CT < 2M when o < and #4,0 > 0. 

Since the NJL model is nonrenormalizable, it is needed to introduce a cutoff in the 
momentum integration. In this paper, we use the three-dimensional momentum cutoff as 

<i 3 P 1 rA 



<*)> ' Jo " (14) 

Hence, the present model has six parameters, m , A, g 2 ,o, g 4j0 , g Qt2 and #2,2- We simply 
assume m = 5.5 MeV. For the case of g^o = 0, A and #2,0 are fixed to reproduce the 
empirical values of the pion decay constant f n = 93.3 MeV and the pion mass M n = 138 
MeV. For the case of nonzero g^ , A, g 2 ,o an d #4,0 are fixed to reproduce f n and M n above 
and M a = 650 MeV. In the latter case, the contribution of the p 4)0 term to G* a is about 



11 percents at T = and fi = 0. For the vector coupling constants #0,2 and #2,2, we take 
two extreme cases, G w = and G u = G a /1.5, where G a = (j* ct |t=^=o = 292,0 + 12(74,000 and 
G^ = G*J T=At=0 = 2g Q 2 + 2(?2,20"o f° r the scalar density at T = and \i = 0. Furthermore, 
in order to determine each value of g ,2 and g 2 ,2, we take two cases, (2g 2 , 2g 2 ,2<3"o) = (G w , 0) 
and (2(yf 0i 2, 2(/2,2C r o) = (0.8(7^, 0.2(2^). For the second case, G* VLU is suppressed in the high \i 
region, as shown later. Table I summarizes the parameter sets we have taken, M a at T = 
and n = 0, and the type of transition at T = 0. 



model 




92,0 


54,0O"0 


2^0,2 


2 fl , 2,2<7o 


M a [GeV] 


type 


NJL 




5.498 











0.681 


first order 


NJL + a 4 




5.276 


0.1109 








0.650 


first order 


NJL + J 2 




5.498 





G<jj 





0.681 


crossover 


NJL + uj 2 


+ a 2 co 2 


5.498 





0.8 G w 


0.2 G, 


0.681 


crossover 


NJL + a 4 


+ u 2 


5.276 


0.1109 


Gu, 





0.650 


crossover 


NJL + a 4 


+ uj 2 + a 2 J 2 


5.276 


0.1109 


0.8 G w 


0.2 G^ 


0.650 


first order 



TABLE I: Summary of the parameter sets, M a at T = and \i = 0, and the type of the 
transition at T = 0. The coupling constants are shown in GeV -2 . The effective coupling G u 
is fixed to 7.331 GeV~ 2 in the NJL+w 2 and the N3L+uj 2 +a 2 u 2 model and to 7.922 GeV" 2 in 
the NJL+<r 4 +w 2 and the NJL+<r 4 +w 2 +cr 2 u; 2 model. For all cases, we take A = 0.6315 GeV and 
M| T=M=0 = 0.3379 GeV and a = -0.03023 GeV 3 . 

Figure 1(a) shows the T dependence of the effective quark mass for the case of fi = 0. 
Since the uo field has no contribution to M when /i = 0, results are shown only for the original 
NJL and the NJL+cx 4 model. The nonlinear a 4 interaction makes the restoration of chiral 
symmetry faster, since the effective coupling, responsible for the SSB of chiral symmetry, 
becomes smaller as T increases, as shown in Fig. 1(b). However, the interaction keeps the 
transition a crossover as in the case of the original NJL. The NJL+cx 4 model yields a smaller 
transition temperature (T c ~180MeV) than the original NJL model does (T c ~190MeV). The 
value, T c ~180MeV, is close to the one predicted by the lattice simulation (T c ~170MeV) 0. 

Figure 2 shows the effective quark mass as a function of fi with T = 0. In the original 
NJL model, the chiral transition is a first order. The nonlinear a 4 interaction makes the 
transition faster, as shown in the result of the NJL+cx 4 model. 
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Fig. 1: The T dependence of (a) the effective quark mass and (b) the effective coupling in the case 
of fi = 0. 
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Fig. 2: The /i dependence of the effective quark mass in the case of T = 0. 
As already shown in Ref. 



16| ]. the uj 2 interaction tends to change the transition from 
a first order to a crossover; actually this is seen in the result of the NJL+co> 2 model of Fig. 
2(a). As an interesting result, however, the nonlinear a 4 and a 2 uj 2 interactions make the 
transition sharper again, so that the transition returns to a first order in the full-fledged 
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NJL+ct 4 +u; 2 +<t 2 u; 2 model. Thus, the nonlinear cr 4 and a 2 uj 2 interactions work so as to cancel 
out the well-known effect of the us 2 interaction. 

Figure 3 shows the effective couplings as functions of \x for the case of T = 0. For all 
models, both G* a and are suppressed in the high /i region, but each model has its own 
\x dependence. In the NJL+a 4 model, G* sa decreases suddenly in the high /i region and then 
approaches the value of 2g 2 ,o, because the cx-dependent part of G* a almost vanishes there. 
Similarly, in the NJL+cr 4 +cj 2 model, G^ a decreases rather rapidly but not suddenly as (i 
increases and approaches the value of 2g 2 ,o- The change in the /i dependence of Gl a from the 
the NJL+<7 4 model to the NJL+<t 4 +u; 2 model comes from the fact that the uj 2 interaction 
makes the phase transition weak. In the NJL+u; 2 +<7 2 u; 2 model, G* a decreases gradually as fi 
increases, because the cu-dependent part of gives a negative contribution to the effective 
coupling. In the full-fledged NJL+cr 4 +u; 2 -|-cr 2 LG> 2 model, G* a is suppressed suddenly by the 
a 4 interaction and then suppressed gradually by the a 2 uj 2 mixing interaction. 

As for G* w shown in Fig. 3(b), one can see a similar sharp suppression but not find any 
gradual decrease. This is understandable from the fact that in Eq. (J10)) the coupling has a 
u-dependent term but no tu-dependent one. As a point to be noted, the sharp suppression 
comes from the a 2 uj 2 interaction for the case of G* , but the a 4 interaction for the case of 
Furthermore, note that the <j 2 lo 2 interaction yields both a gradual decrease of G* a and 
a sharp suppression of G* u . 

The \i dependences of the effective couplings mentioned above are similar to that of the 
chiral symmetry restoration shown in Fig. 2. We can then consider that effects of higher- 
order interactions on the chiral symmetry restoration are described mainly through the 
effective couplings in their \i dependences, although U is also changed by the higher-order 
interactions. 

Figure 4 shows the phase diagram in the fi-T plane. Results are shown for the three 
cases of the original NJL model, the NJL+a 4 model and the NJL+cr 4 +u; 2 -|-cr 2 ti; 2 model in 
which the transitions are first order in the high /x and low T region. First-order transitions 
take place on the curves in Fig. 4, and at the endpoint (/x e , T e ) of each curve the transition 
is changed into a crossover. The NJL+a 4 model yields smaller jj, e and larger T e than the 
original NJL model does; (/x e , T e ) = (308, 54) MeV for the former and (330, 47) MeV for 
the latter. When we take another parameter set, A = 0.6315 GeV, #2,0 = 5.00 GeV -2 and 
£4,0 = 271 GeV" 8 , that reproduces M a = 600 MeV, the endpoint of the NJL+cr 4 model 
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Fig. 3: Effective couplings as functions of /x in the case of T = 0. 




Fig. 4: Phase diagram in the [i-T plane. Each curve denotes the location of the first-order phase 
transition in each model. 

is at (276, 62) MeV. Thus, the nonlinear a interaction shifts the critical endpoint toward 
values predicted by the lattice QCD calculation s, ( u P ,T P ) = (242, 160) MeV [7j, and by the 
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221. 



QCD-like theory, (/x e ,T e ) ~ (200, 100) MeV 

As mentioned above, the uo 2 interaction tends to change the transition from a first order 
to a crossover, but this effect is partially canceled out by the scalar-vector mixing interac- 
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tion a 2 oj 2 . Consequently, as shown in Fig. 4, there exists a critical endpoint also in the 
NJL+<t 4 +w 2 +(t 2 cj 2 model. The comparison of our models with the results of the lattice 
simulations and the QCD-like theory indicates that the NJL+a 4 model is most preferable 
among our models. This may imply that the cancellation is essential. 

In summary, we have studied effects of eight-quark interactions on the chiral phase tran- 
sition. The scalar-type nonlinear interaction a 4 hastens the restoration of chiral symmetry 
and shifts the critical endpoint toward the values predicted by the lattice simulations and 
the QCD-like theory. The vector-scalar mixing interaction a 2 u 2 can make the transition 
sharper in the high /i and low T region, while the u 2 interaction works in the opposite direc- 
tion. Thus, the effect of the mixing interaction tends to cancel out that of the oo 2 interaction 
in the high \x region, while the u 2 interaction is still strong in the normal nuclear density 
region. The roles of the eight-quark interactions are well understood through the effective 
couplings, G* CT , G* CT , G* w and in their /x dependences. Our results indicate that eight- 
quark interactions are very important for phase transitions and must be studied in detail. 
It is very interesting to study effects of the interactions on a color superconductivity itself 
and its correlation with the chiral phase transition. The study is now in progress. 
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